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Abstract 

We present a study of the decay r~ — > Ks~k~v t using a 351 fb _1 data sample 
collected with the Belle detector. The analysis is based on 53, 110 lepton-tagged 
signal events. The measured branching fraction B(t~ — > Kstt~v t ) = (0.404 ± 
0.002(stat.) ±0.013(syst.))% is consistent with the world average value and has bet- 
ter accuracy. An analysis of the Ksir~ invariant mass spectrum reveals contributions 
from the if* (892) ~ as well as other states. For the first time the K* (892) ~ mass 
and width have been measured in r decay: M(K* (892)~) = (895.47 ± 0.20(stat.) ± 
0.44(syst.) ± 0.59(mod.)) MeV/c 2 , r(A"*(892) - ) = (46.2 ± 0.6(stat.) ± 1.0(syst.) ± 
0.7(mod.)) MeV. The if*(892) _ mass is significantly different from the current world 
average value. 
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PACS: 13.30.Eg, 13.35.Dx, 13.66.Jn, 14.40.Ev, 14.60.Fg 
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1 Introduction 



r lepton hadronic decays provide a laboratory for the study of low energy 
hadronic currents under very clean conditions. In these decays, the hadronic 
system is produced from the QCD vacuum via the charged weak current me- 
diated by a W ± boson. The r decay amplitude can thus be factorized into a 
purely leptonic part including the r and v T and a hadronic spectral function. 
Strangeness changing r decays are suppressed by a factor of ~ 20 relative 
to Cabibbo-allowed modes. High-statistics measurements at B factories pro- 
vide excellent opportunities for studying the structure of the strange hadronic 
spectral functions in specific decay modes [1,2,3], the parameters of the inter- 
mediate states and the total strange hadronic spectral function [4]. 

The decay t~ — > K it~v T (unless specified otherwise, charge conjugate decays 
are implied throughout the paper) has the largest branching fraction of all 
Cabibbo-suppressed decays of the r lepton. Early studies of this decay es- 
tablished that the main contribution to the Kn invariant mass spectrum is 
from the fT*(892) meson [5,6,7]. Although scalar or tensor contributions are 
expected in theoretical models [8,9] and not excluded experimentally [10,11], 
the low statistics of previous investigations did not allow for a detailed study. 

Here we report a precise measurement of the branching fraction for the decay 
t~ — > K s it~is t as well as a study of its final state dynamics. This analysis is 
based on a 351 fb _1 data sample that contains 313 xlO 6 t + t _ pairs, collected 
with the Belle detector at the KEKB energy-asymmetric e + e _ (3.5 on 8 GeV) 
collider [12] operating at the T(4S) resonance. 



2 The Belle detector 



The Belle detector is a large-solid-angle magnetic spectrometer that consists 
of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an 
array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrange- 
ment of time-of-flight scintillation counters (TOF), and an electromagnetic 
calorimeter (ECL) comprised of CsI(Tl) crystals located inside a supercon- 
ducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- return 
located outside the coil is instrumented to detect K\ mesons and to identify 
muons (KLM). Two inner detector configurations are used in this analysis. 
A beampipe with a radius of 2.0 cm and a 3-layer silicon vertex detector are 
used for the first sample of 124 xlO 6 t + t~ pairs, while a 1.5 cm beampipe, a 
4-layer silicon detector and a small-cell inner drift chamber are used to record 
the remaining 189 xlO 6 t + t~ pairs [13]. The detector is described in detail 
elsewhere [14]. 
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3 Selection of r + r events 

We select events in which one r decays to leptons, r~ — > l~V\v T , I = e, /i, 
while the other one decays via the hadronic channel t~ — > h~u T , where h~ 
denotes the hadronic system. Events where both r's decay to leptons are used 
for normalization. This reduces systematic uncertainties substantially. 

The selection process, which is optimized to suppress background while re- 
taining a high efficiency for the decays under study, proceeds in two stages. 
The criteria of the first stage suppress beam background to a negligible level 
and reject most of the background from other physical processes. These cri- 
teria retain a 46.0% efficiency for r + r~ events. We then select events having 
2 to 4 tracks with a net charge less than or equal to one in absolute value. 
The extrapolation of each track to the interaction point (IP) is required to 
pass within ±0.5 cm in the transverse direction and ±2.5 cm in the longitu- 
dinal direction of the nominal collision point of the beams. Each track must 
have a transverse momentum in the center-of-mass (CM) frame larger than 
0.1 GeV/c. At least one of the charged particles should have a transverse mo- 
mentum higher than 0.5 GeV/c. The sum of the absolute values of the CM 
track momenta must be less than 9 GeV/c. The minimum opening angle for 
any pair of tracks is required to be larger than 20°. The number of photons 
with a CM energy exceeding 80 MeV is required to be less than or equal to 
five. The total ECL energy deposition in the laboratory frame must be less 
than 9 GeV. The total energy of all photon candidates in the laboratory frame 
should satisfy J^E}j AB < 0.2 GeV. The missing four-momentum P m i SS is cal- 
culated by subtracting the four-momentum of all charged tracks and photons 
from the beam four-momentum. The missing mass M m j SS = \J P^ liss is required 
to satisfy 1 GeV/c 2 < M miss < 7 GeV/c 2 . The polar angle of the missing 
momentum in the CM frame is required to be larger than or equal to 30° and 
less than or equal to 150°. The last two criteria are particularly effective in 
suppressing the backgrounds from radiative Bhabha, e + e~ — > yU + /i~(7) and 
two-photon processes. 

At the second stage, two event classes are selected for further processing: a two- 
lepton sample (if, 1%), l±,l 2 — e, fj, and a lepton-hadron sample (^, K s tc t ), I = 
e, [i. To select electrons, a likelihood ratio requirement V e = C e / (C e +C x ) > 0.8 
is applied, where the electron likelihood function C e and the non-electron func- 
tion C x include information on the specific ionization (dE/dx) measurement 
by the CDC, the ratio of the cluster energy in the ECL to the track momentum 
measured in the CDC, the transverse ECL shower shape and the light yield 
in the ACC [15]. The efficiency of this requirement for electrons is 93%. To 
select muons, a likelihood ratio requirement = C^/ (C^ + C n + Ck) > 0.8 is 
applied. It provides 88% efficiency for muons. Each of the muon(£ M ), pion(£ 7r ) 
and kaon(£x) likelihood functions is evaluated from the information on the 
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difference between the range calculated from the momentum of the particle 
and the range measured by KLM, and the x 2 of the KLM hits with respect to 
the extrapolated track [16]. To separate pions from kaons, for each track we 
determine the pion (C' n ) and kaon (C K ) likelihoods from the ACC response, 
the dE/dx measurement in the CDC and the TOF flight-time measurement, 
and form a likelihood ratio Vx/n = k/{£'-k-\-£' k)- For pions we apply the re- 
quirement Vk/tt < 0.3, which provides a pion identification efficiency of about 
93%, while keeping the pion fake rate at the 6% level. 

To evaluate the background and to calculate efficiencies, a Monte Carlo (MC) 
sample of 1.50 x 10 9 t + t~ pairs is produced with the KORALB/TAUOLA 
generators [17,18]. The detector response is simulated by a GEANT3 based 
program [19]. 

3.1 Two-lepton events 

For this class the (e, e) and (fi, fi) samples still contain contamination from 
radiative Bhabha and e + e~ — > fx + fi~( / -f) processes of about 50%, only (e, //) 
events are used for normalization. To further suppress BB and charm back- 
grounds, we require the opening angle of the leptons to be larger than 90° in 
the CM. As a result, we selected 2,018,000 events of the (e + ,/i~) type and 
2,028,000 (e-,//+) events. 

MC simulation indicates that there is an approximately 5% contamination 
coming primarily from the two-photon process e + e~ — > e + e~fi + fi~ (2.0%) 
and from t + t~ — > e + (fj l + )7r~v e (v ll )u T i7 T events where the n is misidentified as 
a lepton (2.8%). Contamination from other non-r + T~ processes is found to 
be negligible (less than 0.1%). The numbers of (e + ,/i~) and (e~,/i + ) events 
after background subtraction are 1, 929, 300 ± 1, 400 and 1, 911, 700 ± 1, 400, 
respectively. The detection efficiencies and their statistical errors are (19.262 ± 
0.006)% for (e + ,/i-) and (19.252 ± 0.006)% for (e",/i+) events. 

3.2 Lepton-hadron events 

For this class we select events with only one lepton (/ = e, //), one K$ 
candidate and one charged pion 7r T . A K$ meson is reconstructed from a pair 
of oppositely charged pions having invariant mass M nw (K s ) within ±13.5 MeV 
of the K s mass, which corresponds to a ±5a signal range. The pion momenta 
are then refitted with a common vertex constraint. The ^-distance between 
the two helices at the 7r + 7r~ vertex position before the fit is required to be 
less than 1.5 cm, where z is defined as the direction opposite to the positron 
beam. The closest approach of at least one track to the IP in the r — ip plane 
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must be larger than 0.03 cm. The decay length of the Ks candidate in the 
r — tp plane must satisfy 0.1 cm < L± < 20 cm. The ^-projection of the Ks 
candidate decay length is required to be L z < 20 cm. The Ks decay length 
L(Ks) = \J L 2 ± + L\ must be larger than 2 cm. The cosine of the azimuthal 
angle between the momentum vector and the decay vertex vector of the Ks 
candidate is required to be larger than or equal to 0.95. The lepton-i^^ and 
lepton-7r opening angles are required to be larger than 90° in the CM. 68, 107 
events were selected for further analysis. Figure 1 shows a comparison of the 
MC and experimental distributions for the 7r + 7r~ invariant mass of the Ks 
candidate and the K s decay length. 




MJKs), GeV/c 2 




L(Ks), cm 



Fig. 1. MC (histogram) and experimental data (points) distributions normalized 
to the same number of events, (a) shows the 7r + 7r~ invariant mass distribution 
for Ks candidates, (b) shows the Ks candidate decay length. For each distri- 
bution all the criteria described in the text except the one pertaining to the 
displayed parameter are applied. Applied cuts are shown by vertical lines. 



Figure 1 (a) shows that MC 7t + tt~ mass resolution is slightly better than the 
experimental one resulting in a clear difference of the 7r + 7r~ mass spectra in 
the region of the Ks peak. However, the efficiency of the cut for the 
Kg candidates is almost 100%, hence the impact of this discrepancy on the 
detection efficiency is very small and is taken into account in the systematic 
uncertainty. In Fig. 1 (b) one can see a clear difference between the L(Ks) 
distributions in the region of small L(Ks), where events of r~ — > 7r~7r~7r + z/ r 
decay are located, however, in the region, where L(K S ) > 2 cm, populated 
mostly by true Kg's the agreement is good. Figure 2 shows selected events 
on a plot of the Ks decay length versus the 7r + n~ invariant mass of the Ks 
candidate. The main background is from other r decays: r~ — > KsTt~ Klu t , 
t~ — > iCs 71 " 71 " ^) r ~ ~~ * KsK^u T , t~ — > 7t~7t~tt + v t . Using the branching frac- 
tions of these decays from Ref. [20] and detection efficiencies from MC sim- 
ulation, the contamination from decays with a Ks is calculated to be 14.7%. 
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Fig. 2. Decay length vs. ttit invariant mass of the Kg candidate for (e + , KsTt~) 
events. All selection criteria described in the text except for those pertain- 
ing to the parameters being displayed are applied, (a) shows MC data, where 
events with a real K$ are plotted as points, and the events with fake Kg's 
from t~ — > n + ir~ tt~ u T are plotted as boxes, whose sizes are proportional to the 
number of entries, (b) shows experimental data. The signal region is indicated 
by the middle rectangle, while sideband regions are shown by the rectangles to 
the left and right of the signal region. 

t~ — ► 7t~7t~7t + v t decays contaminate the sample when a pair of oppositely 
charged pions is reconstructed as a fake Kg. The tt + tt~ invariant mass dis- 
tribution of these fake Kg's is flat in the region of the Kg mass (see also 
Fig. 2). The number of 3n background events is calculated from two sideband 
regions in the L{Kg) vs. M^{Ks) plane, determined by the following criteria: 
468 MeV/c 2 < M n7T (K s ) < 482 MeV/c 2 and L(K S ) > 2 cm for the first region, 
515 MeV/c 2 < M n7r (K s ) < 528 MeV/c 2 and L(K S ) > 2 cm for the second one. 
These sidebands have the same area as the signal region. The fraction of signal 
events in the 37r-sideband region is about 1%, which is taken into account in 
the calculation of the MC signal detection efficiency. We observe a 5.6% back- 
ground of 3n events in the signal region. In the {l^^Ks^), I = e, fi sample 
there is a small contamination (of about 0.3% for the e-tagged and 2.4% for 
the /i-tagged events) coming primarily from (tt ± , Ks^ t ) events, where the first 
pion was misidentified as a lepton. The non-r + r~ background is found to be 
0.6%. After background subtraction 53, 110±271 signal events remain. Table 1 
shows how they are distributed among the various tagging configurations. 



4 t — > v T branching fraction 



The r — > K$t( v t branching fraction is calculated according to the formula: 
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Table 1 

Branching fractions for different tagging configurations 







(e-,K S Tr + ) 




(fl-,K S TT + ) 


^cxp 


13336 ± 137 


13308 ± 137 


13230 ± 134 


13236 ± 134 


e(l,K s n),% 


5.70 ± 0.02 


5.58 ±0.02 


5.95 ± 0.02 


5.89 ±0.02 


B(K s wu),% 


0.406 ± 0.005 


0.414 ± 0.005 


0.397 ± 0.005 


0.400 ± 0.005 


<B> h % 


0.410 ± 0.003 


0.399 ± 0.003 


< B > a ii, % 


0.404 ± 0.002 



B(K s ^u T ) 



N(t£,Ks**) e(lf,lf) 



B(l^uiu T ), Zi )2 = e,n, (1) 



N(lf,q) e(lf,K s ^) 

where N(lf, K s tt t ), e(lf, K s tc t ) are the number and MC efficiency of the 
signal (if, Kstt t ) events, N(lf, 1%), e(lf, /J) are the number and MC efficiency 
of the two-lepton (if, t^) events, B(lfuiu T ) is the r leptonic branching fraction 
taken from Ref. [20]. Note that the tag-lepton (if) efficiency cancels in the 
ratio of the efficiencies, so the associated systematic uncertainty is reduced. 
The branching fractions calculated separately for each event configuration are 
given in Table 1, which also lists separately the averages for electrons and 
muons as well as the overall branching fraction. 



Table 2 

Systematic uncertainties 


Source 


Contribution,% 


Ks detection efficiency 


2.5 


t + t~ background subtraction 


1.6 




1.0 


Lepton identification efficiency 


0.8 


Pion momentum 


0.5 


Non-r + T~ background subtraction 


0.3 


B(lviv T ) 


0.3 


e(li,Ksir) 


0.2 


Ks momentum 


0.2 


Pion identification efficiency 


0.1 


Total 


3.3 



Table 2 lists the different sources of systematic uncertainties for the branching 
fraction. The dominant contributions come from the Ks detection efficiency 
and background subtraction. A systematic uncertainty in the K s detection 
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efficiency receives contributions from the reconstruction of Kg daughter pions 
(2.3%), the efficiency for fitting two pion tracks to a common 7r + 7r~ vertex 
(0.9%), which was evaluated by varying the cut on the z-distance between the 
two helices at the vertex position before the fit, and the efficiency of the selec- 
tion criteria (0.6%), which was checked by varying cuts on the it + it~ invariant 
mass M n7T (Ks)- Systematic uncertainties arising from T + r~-background sub- 
traction are 0.8%, 1.1%, 0.6% and 0.5% for the r~ -> K s K L n-v T , r~ -> 
K s n~n v ri r~ — > K s K~v T and t~ — > ix~ix~ix + v T modes, respectively. For 
the background from r decay modes with a K s the uncertainties are deter- 
mined by the corresponding uncertainties in their branching fractions taken 
from Ref. [20], except for the t~ — > KsKlh~v t mode. Here we rely on 
the isospin relation B(r~ — > KsKlti^v t ) = 1/2B(t~ — > K + K~n~i' T ) and 
the CLEO result [21] to calculate the t~ — > K s K l t;~v branching fraction 
B(t~ -> K s K l -r-v t ) = (0.078 ± 0.006)%. The uncertainty in the contamina- 
tion by t~ — > 7r _ 7r _ 7r + kv events is evaluated by varying the K s decay length 
cut. 

The lepton detection efficiency is corrected using the e + e~ — > e + e~/ + Z~, / = 
e, /i two-photon data sample. An efficiency correction table is calculated in 70 
bins on the plane of momentum vs. polar angle in the laboratory frame and 
then applied to the Monte Carlo efficiencies e(lf, Ks7i T ) and e(lf,l^). Hence, 
the uncertainty on the leptonic efficiency is determined by the statistics of 
the e + e~ — > e + e~/ + Z~ sample and the long-term stability, which is evaluated 
from the variation of the corrections calculated for time ordered subsamples 
of the experimental two-photon data. The pion identification efficiency in MC 
differs from that in data. In the signal sample, Kg mesons provide a source 
of identified pions, which are used to calculate corrections to the MC effi- 
ciency. Therefore, the systematic uncertainty on the pion identification effi- 
ciency is determined by the statistical error of the correction, which is about 
0.1%. To calculate e(l 1: K s ir) a signal MC sample is produced according to the 
7^(892) + fT*(1680) model and the model dependence of e(h, K s n) is found 
to be negligible. 

We also vary cuts on the pion momentum, the kaon momentum, and the total 
laboratory energy of photons {J2 E^ AB ) to check the stability of the branching 
fraction. The total systematic uncertainty of 3.3% is obtained by adding all 
the contributions in quadrature. Our final result for the branching fraction is 
B(t~ -> K s n-u T ) = (0.404 ± 0.002(stat.) ± 0.013(syst.))%. 



5 Analysis of the r — > K S 7r v T spectrum 

The Ks7T~ invariant mass distribution shown in Fig. 3 exhibits a very clear 
if* (892) ~ signal. We parameterize this spectrum by the following function 
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(see Ref. [8] for more detail): 



where s is the K s n~ invariant mass squared and P is the K s momentum in 
the K s ti~ rest frame: 



P(s) 



2v^ 



s — [m K + m w ) 2 s — [m K — m n ) 



(3) 



The vector form factor F v is parameterized by the if* (892), if* (1410) and 
if* (1680) meson amplitudes: 



V 



l + P + X 



BW K *(892){s) + /3BW K *( U10) (s) + xBW K * (1680) (s) 



, (4) 



where f3 and x are complex coefficients for the fractions of the if* (1410) and 
if*(1680) resonances, respectively. BW R (s), (R = If* (892), if*(1410), if*(1680)) 
is a relativistic Breit-Wigner function: 



M fl (s) 



Ml 



s-M 2 R + i^7sT R ( S y 
where Fr(s) is the s-dependent total width of the resonance: 



(5) 



r fl (s) = r, 



M 



ii 



g(g) 

P(M 2 R ) 



21+1 



(6) 



where £ = 1(0) if the Kn system originates in the P(S')-wave state and T 0R is 
the resonance width at its peak. 

The scalar form factor F$ includes the if,* (800) and 1^(1430) contributions, 
their fractions are described respectively by the complex constants h and 7: 



F s = x- 



s s 
— 2 BW Ksm) (s) + 7—2 BW KS{M30) (s). (7) 

1VI K*(800) iW X*(1430) 



The experimental distribution is approximated in the mass range from 0.63 GeV/c 2 
to 1.78 GeV/c 2 by a function calculated from the convolution of the spectrum 
given by Eq. (2) and the detector response function, which takes into account 
the efficiency and finite resolution of the detector. In all fits the if* (892) mass 
and width as well as the total normalization are free parameters. Only the 
strengths (fractions) of the other if*'s are free parameters, while their masses 
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and widths are fixed at the world average values [20]. In the approximation 
x is chosen to be real, because Fg is defined up to the common phase, which 
cancels in \Fg\ 2 . 




Fig. 3. Comparison of the Kgir mass distributions, points are experimental 
data, histograms are spectra expected for different models, (a) shows the fitted 
result with the model incorporating the IT* (892) alone, here the background has 
been already subtracted from both experimental and expected distributions, 
(b) shows the fitted result with the K* (892) + iC * (800) + K*( 1410) model, here 
different types of background are also shown. 

Figure 3 (a) and Table 3 show that the if* (892) alone is not sufficient to 
describe the Kgir mass spectrum. To describe the enhancement near threshold, 
we introduce a Kg (800) amplitude, while for description of the distribution 
at higher invariant masses we try to include the if*(1410), if* (1680) vector 
resonances (see Table 3) or the scalar ^(1430) (see Table 4). Figure 3 (b) 
demonstrates the good quality of the fit with the if<*(800)+if*(892)+if*(1410) 
model. It can be seen from Tables 3, 4 that we cannot distinguish between the 
if*(800) + if*(892) + if*(1410) and #5(800) + if*(892) + if*(1430) models. 
The fit quality with the if *(800) + if* (892) + if* (1680) model (see the fourth 
column of Table 3) is worse than that of the if *(800) + if* (892) + if* (1410) 
and if *(800) + if* (892) + if* (1430) models. 

It should be noted that the absolute value of a sum of two Breit-Wigner 
functions of mass (y/s) can have the same shape for two different sets of 
parameters. In the case of the if *(800) + if* (892) + if* (1430) model the 
relevant parameters are x, I7I and arg(7). This statement holds true when 
mass-independent widths are considered. If the width is mass- dependent, some 
difference in the spectra appears. If in the fit to the data the errors are 
large enough, we cannot distinguish these solutions by their x 2 values. For 
high statistics the two solutions can be distinguished by a x 2 test- While 
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Table 3 

Results of the fit of the Ksn mass spectrum in different models of the non-K *(892) 
mechanism: the IT* (1410) and K*(1680) contributions are described by the complex 
constants (3 and %, respectively, while that from the -fTg(800) is described by the 
real constant x. Masses and widths of the non-.fr* (892) resonances were fixed at 
their PDG values (the Kq(800) ma ss and width were fixed from Ref. [22]). 









KZ(ROQ) 4- /<'*('892">-l- 
+if*(1410) 


KX(R00) -4- K*(RQ1)-i- 
+i^*(1680) 




MpV /r 2 


QQK tco l A in 

Ot?U.<JO _l_ 


one 17 x n on 


004 00 4. n on 


it* (892) i 


MeV 


49.29 ± 0.46 


46.19 ± 0.57 


45.52 ± 0.51 


1/31 






075 ± 006 




arg(/3) 






1.44 ±0.15 




Ixl 








0.017 

0.117 ± 

0.033 


arg(x) 








3.17 ±0.47 








1.57 ±0.23 


1.53 ±0.24 


X 2 /n.d.f. 




448.4/87 


90.2/84 


106.8/84 


P(x 2 ),% 







30 


5 



Table 4 

Results of the fit of the K s tt mass spectrum in the # *(800) + K*(892) + #5(1430) 
model (two solutions). The Xq(1430) contribution is described by the complex con- 
stant 7, while that from the Kq(800) is described by the real constant x. Masses and 
widths of the non-i^*(892) resonances were fixed at their PDG values (the if,* (800) 
mass and width were fixed from Ref. [22]). 

lfo(800) + K*(892) + #*(1430) 





solution 1 


solution 2 


M K*(892)- 


, MeV/c 2 


895.42 ±0.19 


895.50 ±0.22 


r_ftT*(892)- > 


MeV 


46.14 ±0.55 


46.20 ±0.69 


l7l 




0.954 ±0.081 


1.92 ±0.20 


arg( 7 ) 




0.62 ±0.34 


4.03 ± 0.09 






1.27 ±0.22 


2.28 ± 0.47 


X 2 /n.d.f. 




86.5/84 


95.1/84 


P(x 2 ),% 




41 


19 



for the #3(800) + #*(892) + #*(1410) and #*(800) + #*(892) + #*(1680) 
models with a complicated vector form factor the values of \ 2 are signifi- 
cantly different (due to the small (~ 1%) errors at the #*(892) peak), in 
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Table 5 

Results of the fit of the Kstr mass spectrum in the model when the non-K*(892) 
mechanism is introduced by the LASS scalar form factor, described by the param- 
eters a and b. 







K*(892)+LASS 


K*(892)±LASS 






a, b - fixed 


a, b - free 


M K*(892)-i 


MeV/c 2 


895.42 ±0.19 


895.38 ±0.23 


r i^*(892)-) 


MeV 


46.46 ± 0.47 


46.53 ± 0.50 


A 




0.282 ±0.011 


0.298 ±0.012 


a, (GeV/c) 


-l 


2.13 ±0.10 


7.4 

10.9 ± 

3.0 


b, (GeV/c) 


-l 


3.96 ±0.31 


4.5 

19.0 ± 

3.6 


rVn.d.f. 




196.9/86 


97.3/83 


P(x 2 ),% 




io- 8 
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the i^ *(800) + if* (892) + 7^(1430) case with a complicated scalar form fac- 
tor different solutions result in similar P(x 2 ) values (see Table 4) due to the 
relatively low statistics in the region of the Kq(800) and 1^(1430) peaks. 

An alternative way to describe our data is to use the parameterization of the 
scalar form factor suggested by the LASS experiment [23,24]: 

Fs = \A LASS (s), A LASS = ^-(sm5 B e iSB + e 2lS * BW K * {im (s)), (8) 

where A is a real constant, P is K$ momentum in the Ksir rest frame (see 
Eq. (3)), and the phase 8b is determined from the equation cotJg = ^ + 
where a, b are the model parameters. In this parameterization the non-resonant 
mechanism is given by the effective range term sin5se l<5s , while the resonant 
structure is described by the 1^(1430) amplitude. 

Table 5 shows the results of fits to the spectrum in models, where the non- 
-fT*(892) mechanism is described by the LASS parameterization of the scalar 
form factor. In the first fit (see the second column of Table 5) a and b pa- 
rameters were fixed at the LASS optimal values [24]. In the second fit a and 
b were free parameters (see the third column of Table 5). The optimal values 
of a and b in our fit differ significantly from the values obtained by the LASS 
collaboration in experiments on Kn scattering [23]. 

The K*{800) + K*{892) + if* (1410) model was considered as the default and 
was used to obtain the K* (892) (K s rc)iy fraction in the K s -kv final state, which 
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was found to be B(t~ -> K\892)~v T ) ■ B(K*(892)- -> K s tt-)/B(t- -> 
Ks7T~u T ) = 0.933 ± 0.027. The 0.027 error includes the model uncertainty, 
which was found by calculating this fraction in the fits with the other models 
mentioned above, as well as the uncertainty in the fit parameters. Finally we 
obtain B(t~ -> K*(892)-u T ) ■ B(K*(892)~ -> AVO = (3.77 ± 0.02(stat.) ± 
0.12(syst.) ± 0.12(mod.)) x 10~ 3 . 



6 Measurement of the A'* (892) mass and width 

A fit to the K s n~ invariant mass spectrum also provides a high precision 
measurement of the A'*(892)~ mass and width. We consider a fit with the 
^(800) + A^*(892) + A"*(1410) model, which provides a good description of 
the data, as a reference, and use it to obtain the A"*(892) - mass and width 
values. It can be seen from Table 3 that the statistical uncertainty is about 
0.20 MeV/c 2 for the mass and 0.6 MeV for the width. Two additional sources 
of uncertainty are studied: the effects of imperfect knowledge of the detector 
response function and model uncertainty. 

The systematic uncertainty is studied with a MC sample by comparing the 
A'*(892)~ parameters implemented in the generator and its parameters after 
the full reconstruction procedure (the detector response function is determined 
from other statistically independent MC simulations of signal events). It is 
found to be 0.44 MeV/c 2 for the mass and 1.0 MeV for the width. 

The model uncertainty is investigated by fitting the K s ti~ mass spectrum with 
different models. The maximal difference from the reference value is considered 
as a model uncertainty. It is found to be 0.59 MeV/c 2 for the mass and 0.7 MeV 
for the width. 

As a result, the AT*(892)~ mass and width are M(AT*(892)~) = (895.47 ± 
0.20(stat.) ± 0.44(syst.) ± 0.59(mod.)) MeV/c 2 and r(AT*(892)-) = (46.2 ± 
0.6(stat.) ± 1.0(syst.) ±0.7(mod.)) MeV, where the first uncertainty is statis- 
tical, the second is systematic and the third is from the model. 



7 Conclusions 

The branching fraction of the r~ — > Ksts~v t decay has been measured using 
a data sample of 351.4 fb _1 collected with the Belle detector. Our result is: 

B{j- -> K s ix-v T ) = (0.404 ± 0.002(stat.) ± 0.013(syst.))% 



15 



To compare our result with the previous measurements made by the OPAL [27], 
ALEPH [10,28], CLEO [11] and L3 [29] groups we calculate the t~ -> K°n-v T 
branching fraction according to the formula B(t~ — > K tt~v t ) = B(t~ — > 
Ks7T~u T ) + B{t~ — > K L it~v T ) = 2B{t~ — > K s it~v T ) and obtain: 

B{j- -> K*ix-v T ) = (0.808 ± 0.004(stat.) ± 0.026(syst.))% 

Figure 4 (a) shows the results of various measurements of the r~ — > K°7T~v T 
branching fraction, along with the Particle Data Group (PDG) fit value (Bpdg( t 
K°tt-v t ) = (0.900 ±0.040)%) [20] and our result. Our result is consistent with 
previous measurements, but is more precise. 



a) 



Belle 



PDG07FIT 



OPAL-00 



ALEPH-99 
ALEPH-98 



1 

CLEO-96 



L3-95 



i i i i I i i i i 



0.6 0.7 0.8 „0.9 



Br(K°nv),% 



1.1 



K (892) 



PDG07 I 

K"%892) 

\-*-\Beltie 




885 887.5 890 892.5 895 897.5 

m k- ( 892)> MeV/c 



Fig. 4. Comparison of the t~ — > K°ir~v T branching fraction (a) and X*(892) _ 
mass (b) measured in different experiments, (b) also shows all available data on 
the X*(892)~ mass together with the PDG average (the hatched region marks 
PDG data, which were not used in the calculation of the average mass, see 
Ref. [20]), as well as our result, which is close to the PDG i^*(892)° mass. 

The i^*(892) alone is not sufficient to describe the Kstt invariant mass spec- 
trum. The best description is achieved in the #5(800) + #*(892) + #5(1410) 
and # *(800) + #*(892) + #5(1430) models. Future high precision studies 
of the invariant mass spectra in r lepton decays with kaons combined with 
angular analysis, i.e. an application of the structure function formalism sug- 
gested in Ref. [1], will elucidate the nature of the scalar form factor. They 
will also check various theoretical models describing the scalar #7r sector, 
e.g., the predictions of the resonance chiral theory [25] and the parameters of 
the #5(800) resonance calculated from the Roy-Steiner representations in a 
model- independent way [26]. 

The product of r~ — > #*(892)~i/ r and #*(892) _ — > K s n~ branching fractions 
is found to be: 
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B(t~ -> K*{892)~u T ) ■ B(K*(892)- -> i^Tr") = 
(3.77 ± 0.02(stat.) ± 0.12(syst.) ± 0.12(mod.)) x lO" 3 , 

also the if* (892) ~ mass and width are measured: 

M(K*(892)~) = (895.47 ± 0.20(stat.) ±0.44(syst.) ±0.59(mod.)) MeV/c 2 

r(X*(892)-) = (46.2 ± 0.6(stat.) ± 1.0(syst.) ± 0.7(mod.)) MeV 

The values of the if* (892) ~ mass and width that we obtain are more precise 
than any of the existing measurements of these quantities listed in Ref. [20] 
and shown in Fig. 4 (b). While our determination of the width is compatible 
with most of the previous measurements within experimental errors, our mass 
value is systematically higher than those before and is in fact consistent with 
the world average value of the neutral if* (892)° mass, which is (896.00 ± 
0.25) MeV/c 2 [20]. Note that all earlier mass measurements listed in Ref. [20] 
come from analysis of hadronic reactions and include the effects of final state 
interaction while our work presents a measurement based on t~ decays, where 
the decay products of the if* (892) _ are the only hadrons involved. It is also 
noteworthy that none of the previous measurements in Ref. [20], all of which 
were performed more than 20 years ago, present the systematic uncertainties 
for their measurements. Unfortunately, previous studies of the if* (892) ~ in 
t~ lepton decays usually do not determine its parameters. The only published 
result we are aware of is that of ALEPH [30], which is consistent with ours. 
Its accuracy, however, is much worse and no systematic errors are presented, 
which precludes any detailed comparisons. A similar if* (892) ~ mass shift of 
(+4.7 ± 0.9) MeV/c 2 was reported by CLEO [31], but no dedicated study of 
this effect was published. Future dedicated measurements of the if* (892) ~ 
parameters with high precision are necessary to clarify this discrepancy and 
shed light on the long standing issue of the electromagnetic mass difference 
between the charged and neutral if* (892) [32,33]. 
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